Heun Functions and the energy spectrum of a 
charged particle on a sphere under magnetic field and 

Coulomb force 

A.Ralko and T.T.Truong, 
Laboratoire de Physique Theorique et Modelisation, 
Universite de Cergy-Pontoise, F-95031, Cergy-Pontoise Cedex, France 

February 1, 2008 

Abstract 

We study the competitive action of magnetic field, Coulomb repulsion and 
space curvature on the motion of a charged particle. The three types of interaction 
are characterized by three basic lengths: Ib the magnetic length, Iq the Bohr 
radius and R the radius of the sphere. The energy spectrum of the particle is 
found by solving a Schrodinger equation of the Heun type, using the technique 
of continued fractions. It displays a rich set of functioning regimes where ratios 
■p- and p take definite values. 
Pacs: 03.65.Ge, 73.43.-f 



1 Introduction 

The motion of charged particles under constant magnetic field on a sphere has been 
introduced by F.M.Haldane [Q and several authors H, Q to discuss the fractional 
quantum Hall effect, and compute several relevant quantities for a system of N inde- 
pendent particles. The sphere is basically a macroscopic support for the motion of a 
large assembly of particles and its radius R is usually of macroscopic order. 

However in recent years, it has been possible to realize spheres of nanostructures 
ranging from fuUerenes size to spherical nanosize objects such as Si02 balls in Opals 

p|. The motion of a charged particle on such spheres with constant magnetic field 
parallel to a given diameter of the sphere has been studied quite extensively [@, ^ . 
The energy spectrum helps to understand some interesting aspects such as orbital 
magnetism. 

In view of future developments, it appears useful to consider the motion of charged 
particles on a sphere under the action of competiting forces. Besides the curvature of 
the configuration space and the confining effect of the radial magnetic field, it may be 
interesting to introduce the Coulomb electric force due to an ion positioned somewhere 
on the sphere. This situation may arise in the process of fabricating the actual sphere 
when a charged impurity may slip into the spherical surface and become a center of 
interaction JTo[ . The point is thus to see how such a charged impurity affects the energy 
spectrum whether or not such effects are in practice relevant. Moreover, the choice of 
a magnetic monopole generating a radial magnetic field is made to obtain analytical 
solutions. Since a constant magnetic field parallel to a given diameter does not allow 
us to find an exact wave function, we replace it with a radial field which competing 
effects are well described. The strength and range of these effects can be described by 
three characteristic lengths: 
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Ib — \/ • magnetic confinement length. 



• /o = Y jJkc^ • Bohr radius of the hydrogen atom. 

• R: radius of the sphere. 

Here M and e are respectively the mass and the charge of the particle, B the 
constant radial magnetic field on the sphere and k = -i^— in MKSA unit. 

Problems with three characteristic lengths do occur in physics. We may cite the 
case of one dimensional superconductivity in the presence of twinning defects where 
electron pairs are moving in constant magnetic field and a finite Dirac Comb 
potential. The three characteristic lengths involved are: 

• Ic- coherence length of the electron pairs, 
magnetic confining length. 

• d: period of the Dirac Comb potential. 

In the present case we look for stationary states of a charged particle on a sphere 
under a constant radial magnetic field and repelled by a point charge placed on the 
sphere. This is done by solving the time independent Schrodinger equation. The so- 
lution will turn out to be expressible in terms of a generalization of Gauss Hyperge- 
ometric Function called the Hcun Function This new function has been studied 
by several authors ^ |l4| and only one type of Heun function turns out to be 
relevant for our problem. In fact to fulfill the requirements of quantum mechanics it 
is necessary to demand that the Heun function be defined all over the sphere and also 
of square integrable class. This yields a condition on one of the coefficient of the Hcun 
equation, which leads to the quantization of the energy levels. 

In section 2 we present the classical aspects of the dynamics showing that in partic- 
ular the motion in the ^-direction the energy should be larger than a certain limit. In 
section 3 we introduce the quantum treatment which leads to the Heun equation. We 
compute in section 4 the i? — > c» and S oo limit, which describes the recovery of 
the motion in planar geometry. In this case, the limiting problem is that of the relative 
motion of two equally charged particles on plane in uniform perpendicular magnetic 
field and under Coulomb repulsion. Results and comments are left to section 5. 



2 Classical considerations 

Let {9, (j)) be the angular coordinates of a particle of mass M and charge e on a sphere 
of radius R. There exists a constant radial magnetic field B pointing outward on the 
surface of the sphere. B is described by the vector potential 

Q ^ 

^ = -$0 cot 610 (1) 

where S (the half magnetic flux) is the strength of the magnetic monopole generating 
the magnetic field B and $0 = 7, the elementary flux unit, and cf) the unit vector along 
the (j) direction on the sphere. In fact we have 

V X A = B = [Br = B,Be = B^ = Q). (2) 
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We assume in addition the existence of a repulsive Coulomb interaction due to a 
charge e, placed at the North pole of the sphere, which is given by the potential 

V{9)=n-/^. (3) 
2R sm I 

The dynamics of the particle is given by the Hamiltonian: 

^ pI , ip<p + hScose)^ 

2MR^~^ 2MR^sinH 2i?sin|' 

where pg and are the canonical conjugate momenta to 9 and (j). Motion at constant 
energy is defined by 

H = E. (5) 

For the sake of simplicity of the discussion we shall work with dimcnsionlcss quan- 
tities: He, n^, and e defined by 



Introducing the Bohr radius 



the equations of energy conservation becomes: 



, {U^ + ScosB)^ R 1 

lie H TT^ ^ 1 — — - e- (») 

sm 9 iQ sm 
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We see the competing roles of geometry and Coulomb repulsion expressed by the 
ratio ^. 

Moreover in spherical coordinates we have: 

Br = - \B\ = B. (9) 

we can then introduce a magnetic length Ib = \J and express S as: 

Hence the combined interactions can be expressed through the ratios of lengths 
^ = \/^and f . 

IB to 

Note that since < 9 < n, sin | is always positive and the left hand side of eq.(||) 
is the sum of three positive terms. 

Now since p^ is a conserved quantity (0 is a cyclic variable) one may set p^ = mh, 
with m fixed equal to H^. The orbit in the phase space (Hg, 9) is a curve of equation: 

y R I 



^9 + —^[m+-f^cos9] +-—^=e. (11) 
sm^ 9 \ Ib J k sm T 
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Equation (^Jj) can be viewed as the motion of a fictitious one dimensional particle in 
0-space. The potential function for < 9 < tt has a minimum value cq which depends 
on m, and j^, thus the motion is only possible if e > eg. 
The action for a complete cycle is Jg given by: 



27r 



1 / i?2 \^ i? 1 







Je = 2h dOJe- —^[m+-^cose] - (12) 
' y sin^ 6* V ^1 / ^0 sin f 



Bohr Sommerfeld quantization of the 0-motion by setting: 



Je==/i(7i+i) 71 eN (13) 



yields the quantized energy levels through e = En^m such that: 



Ui+ COS0 _ d0^7r(n+^). (14) 

sm V 's / 'o sm | i 



Unfortunately the integral can not be evaluated in closed form in terms of simple 
functions. This aspect will not be treated here. 



3 Quantum motion 



The Hamiltonian of the motion can be suggestively reexpressed in terms of the oper- 
ators A — rx (p + eA) as in [|l|: 



H = 



A2 



2A/i?2 2i?sinJ 



(15) 



A has cartesian components: 

Ay 

A, 



+ hS cos (j)^: 

My + hS sin (jjf. 
A-f, - hScost 



(16) 



with usual definition of Af : 

My 



ih 



sin (f)-^ + cos (j) cot 9-^ 

— cos (f)-^ + sin (j) cot 
d 



(17) 



consequently with 



.2 I 



1 



one has: 



A^ = Af 2 



52 + 2^5* cose. 



2 c2 



^5 



(18) 



(19) 



and the Schrodinger equation for stationary states ^^[9, (j)) reads: 

1^ - S'2 - 2^5* cose. 



^+cote|)*(e,, 



sm 
R 1 



^0 si 



sm -, 



*( 



(20) 
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As (/) is a cyclic variable, is conserved, thus wave function '^'{9, (f)) can be set of 
the form : 



(21) 



with m e Z and we get a new equation for F{9): 



S^~2mScos9 R 1 



sm 



to sm ; 



F"{9)+cot 9F'{9) 
+ e + 54 F{9) = 0. 



(22) 



The choice of m as a relative integer guaranties the uniformity of the wave function 
under (/(-rotation. We seek solutions F{9) which vanishes at 6* — > because of the 
Coulomb repulsion at the Norh pole and which is square integrable. 

We now transform this equation with the change of variable z — cos 9 to bring it 
back to a known canonical form. F is now a function of z, satisfying: 



+ 



52 - 2mSz V2R 1 



1-zJ 



Iq \J\ — z 



{1- z^)F"{z)-2zF'{z) 
+ (e + S^)\F{z)^0. 



A change of unknown function through the substitution: 

F{z) = (l-z)t(l + z)5F(z), 



(23) 



(24) 



where a and b are free parameters to be chosen later, and a last change of variable 
X = \J^^ transforms equation ( p3| ) into a Heun equation under its canonical form: 

{x - l){x ~ a')P" [x) 



{t' + 6 + 7)x - (e' + a' 6 + 7(1 + a')) + ^ 



P'{x) 



a(3- 



a(3h 



P{x) = 0, 



(25) 



with the relation a + (3 + \ = 7 + (5 + e'. 

This equation is a second-order equation with four regular singularities, located 
at (x= 0, 1, a' and 00). h is called "auxilliary parameter" and is significant for our 
problem. The solution to this Hcun equation is a local Heun function admitting an 
expansion around .t = 0. Such a solution has a convergent expansion in a domain 
centered at a: = and for \x\ < 1 and |a;| < \a'\. The South pole which corresponds to 
X — 1 will be a singular point since the expansion is divergent there. 

To get eg. (p5|) , the coefficients a and b must fulfill the conditions: 



-ASm 



(26) 



As F{9) must not have singularity as 6* ^ (or 9 n) and equivalently P{x) must 
be regular for x — s- (or x — > 1), one must choose a and b positive, i.e. with S* = ■§-: 



a = \S — m\ 



\S- 



(27) 
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The parameters of the Heun equation are then related to the parameters of the 
problem through: 



7 
S 



-1 

2a + 1 
e' 6 + 1 



a/3 = (a + 6)(a + 6 + 2)-4(e + S'2) 
R 

aPh = -4- 



(28) 



One possibility to make our solution valid everywhere on the sphere is to reduce the 
Heun functions down to a Polynomial. But this procedure imposes a condition which 
makes the radius of the sphere discrete. This is, in general physically not acceptable. 

The other possibility is to turn to the procedure of augmented convergence ||l2| 
which makes use an expansion in terms of Hypergeometric functions | p^ , [l^ . Following 
Ronveaux, we first expand our local Heun function P{x) in terms of hypergeometric 
functions 



(29) 



where F[—v, + 5 + 7 — 1; 7; x) is a local solution of a hypergeometric equation which 
matches Heun's solution at the singularities and 1 (same exponents at the singulari- 
ties and 1). As P{x) obeys the Heun equation, the fulfill the three-way recursion 
relation Ilia: 



where 



{v + a- l)(i/-|-7- l)(t/ + w- 1) 



(30) 



(2j/ + u; - l)(2i/ + w - 2) 

Lu — af3h + a v{v + w) 
e'v{v + u))(7 - (5) + [v{v + w) + ap\ [2v{v + w) + j{w - 1)] 
{2iy + w- l){2iy + w + 1) 
{i' + + w-a + + w- (3 + + S) 



{2iy + w + l){2u + w + 2) 



(31) 



with w = 7 + (5 — 1. 

The behavior of {y^} is given by |13|: 













l-X 



(32) 



where X — \/l — x~'^ . 

To determine this limit, we must find the asymptotic behavior of the {c^}. For 
V — > 00, the coefficients of the recursion relation become: 



1 



,1 



(l^-a'K 



1 



(33) 



and eq.(pO|) tends to the critical equation + 2 (2a' — l)p + 1 = (where c^+i has 
been replaced by p*) which admits the roots pi and p2 given by: 



I- A 



\Pi\ 



\P2\ = 



\Pl\ 



(34) 
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where A — vl 



Erdelyi pM has shown that the convergence of the series occurs 



in a region determined by 



X 



1-X 



< 



\pl\ 



\P2\ 



(35) 



This is the so-called augmented convergence phenomena. This inequality implies that 
the point x lies in the interior of the ellipse £ with foci at and 1. Thus our wave 
function P{x) is now well defined on each pole of the sphere. 

The condition for which the solution of the Heun equation has augmented conver- 
gence at a; = and x = 1 (poles of the sphere) is that the parameter h must be a zero 
of the continuous fraction obtained from the three way recursion relation [O, i.e.: 



So 



where 



A,. 



Bi 



= 0, 



B2 + 



B3 + . 



and B,, 



(36) 



(37) 



Hence for h — hn 'we can construct a bona-fide wave function for a stationary state. 
The set of /i„ can be obtained by numerical treatment of the continuous fraction. We 
use the so-called modified Lentz's method [|l6[, calculating by iterations the values of 
the continuous fraction and stopping when a given precision is obtained. 
The energy spectrum of the problem is of the form: 



Iq hn 



a + b 



1 - 



'■B 



(38) 



with n ~ l,2,...cx3. Now for S = jr- and m given, all the parameters of the Heun 

^ . . . 

equation are defined and we can find, with classical numerical methods, the solution 

in the form of the transcendental eq.(|3^). Note that since we are interested by the 
energy spectrum, we have replaced hn in eq.(pq) by the following expression: 



h„ = 



4i? 



k 4.{en + ^)-ia + b)ia + b + 2) 



(39) 



which allows to compute directly the energy values e„. This is done in section 5. 



4 The R,S ^cx) limit 

It is interesting to look at the R,S-^oo limit, since it is relevant for many problems 
in physics. In the fractional quantum Hall effect, a lot of calculations were made in 
spherical ( for simplicity) geometry |l[ then the physical effects studied in planar 
geometry. For example, the computation of the energy of the neutralizing background, 
takes a very simple form in the spherical geometry, since it amounts to put a particle 
of opposite charge at the center of the sphere . In this limit, the curvature of the sphere 
vanishes and we expect to recover the plane. We know that S and R are related by 
equation ([To|), thus when i? — > oo, 5 must go also to infinity. For physical reasons, 

5 must be infinite, because there is a constant perpendicular magnetic field ( i.e. a 
non-vanishing flux) on the infinite plane. 
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Starting with the Schrodinger equation on the sphere, after the substitution (p^), 
one has: 



(1 - z^)P"{z) + [b-a-{a + b + 2)z] P'{z) 



a + b\ / a + b 



1 - 



V2R 



loy/l - Z 



P{z) = 0. 



(40) 



Thus we look for a change of variable to describe the plane. Such a variable can be 
X — 2S'(1 — z) and transforms equation (Eol) into: 



-45* ■ 



a + b 
AS 

R 



x) P'{x) 



-{xP"ix) + ^xP'{x)}^0- 



(41) 



For 5 — > oo, we can neglected the last term of this equation and obtain the asymp- 
totic form of the equation: 



XPasix) 



a + 1 



45 



PUx) 



e + (2±^) (2i± + l) R 



45 



2V5VX 



Pas{x)^0- 



(42) 



With the last change of variable — \[\ ^'^id the fact that for 5 ^ 00, a = 
|5 — m| — > 5, 6 = |5 + m| — > 5 , equation (B3) reduces to a Biconfluent Heun equation 
in his canonical form: 



(5 + /?(l + a)" 



(7 -a- 2)^ 



p<.s{i) = 0, 



(43) 



where a, 7 and 8 are expressed in terms of e, 5 and ^. As shown in [Q, this 
equation describes the relative motion of two equal planar charged particles in uniform 
perpendicular magnetic field and under the Coulomb repulsion. So the link between 
problem and planar problem is established and is locally acceptable when i?, 5 — > 00. 



5 Results and Comments 

We present the results by giving the two first energy levels n = 1, 2 as a function of 
S which is basically proportional to the constant magnetic field B. The value of the 
ratio i?//o is fixed for a given radius R — lOOZg . Figures 1 and 2 give the behavior of 
the ground state energy ei.m as function of the applied magnetic field for positive and 
negative values of m. 

We see clearly that these values are always bounded below as seen from classical 
considerations. For negative m values the curves exhibit a sharp discontinuity in the 
slope of the tangent at some values of the magnetic field (|to| = 5), there is also a 
milder discontinuity in the derivative without change of sign for positive values of m. 

In figure 3, the two situations are put together to display the phenomena of level 
crossing with respect to the magnetic field. 
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Figure 1: Energy for n = 1 versus S (proportional to B) for positive values of m 



The behavior of the level £2,711 for positive values of m is given in figure 4. 

Finally crossovers between levels ei and £2 for positive m are displayed in figure 
5. Note that at i? = we have the spectrum of the charged particle in the presence 
of pure repulsion centered at the North pole. For strong magnetic field, B 00, we 
can neglected the Coulomb repulsion, thus the levels tends to the Landau levels as 
expected. 

Higher levels {n > 2) exhibit similar behavior for m > and to < but calculations 
are more involved so they will not be displayed here. 

Let us examine now the limit of vanishing Coulomb repulsion (/q 00). The energy 
levels of eq.(p8|) become 



1 - 



(44) 



and are no longer given by the They arc to be compared with the Landau levels 
on the sphere given by 



a + b 



(45) 



where n is the Landau level index. We recover then the Lowest Landau level (n = 0) 
on the sphere as exected before. It is important to note that Haldane has considered 
only states with a ^ S — m, b = S + m and < m < S . 



To our knowledge, this is the first instance where a special class of Heun function, 
solution of a linear differential equation with four singularities occurs in quantum 
mechanics. This class of function can handle physical problems which carry three 
characteristic lengths. The spectrum energy presents different regimes, going from the 
pure Coulomb problem (B = 0) to the Landau problem {B 00), passing through 
the well known level mixing regime. It is the main difference with the pure magnetic 
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o — o m= 


1 


□ — □ m= 


2 


— ^ m= 


8 


A — A m= 


9 
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Figure 2: Energy for n = 1 versus S (proportional to B) for negative values of m 



problem 0, |[ ^ where, for weak fields, the energy levels are equally separated (Landau 
levels) instead of the Coulombic behavior encountered here. Thus the impurity plays an 
important role in the weak field regime but not in the high field regime. So, the impurity 
can be ignored for strong magnetic fields, i.e. for S beyond the sharp discontinuity at 
fixed m. The complete study of this wave function and its properties are left for a 
future work. 
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Figure 4: Energy for n = 2 versus S (proportional to B) for positive values of m 
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